Wigner-Mott scaling of transport near 
the two-dimensional metal-insulator transition 
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Electron-electron scattering usually dominates the transport in strongly correlated materials. 
It typically leads to pronounced resistivity maxima in the incoherent regime around the coherence 
temperature T* , reflecting the tendency of carriers to undergo Mott localization following the demise 
of the Fermi liquid. This behavior is best pronounced in the vicinity of interaction-driven (Mott-like) 
metal-insulator transitions, where the T* decreases, while the resistivity maximum p max increases. 
Here we show that, in this regime, the entire family of resistivity curves displays a characteristic 
scaling behavior p(T)/p max fts F(T/T max ), while the p max and T max ~ T* assume a powerlaw 
dependence on the quasi-particle effective mass m*. Remarkably, precisely such trends are found 
from an appropriate scaling analysis of experimental data obtained from diluted two-dimensional 
electron gases in zero magnetic fields. Our analysis provides strong evidence that inelastic electron- 
electron scattering - and not disorder effects - dominates finite temperature transport in these 
systems, validating the Wigner-Mott picture of the two-dimensional metal-insulator transition. 

PACS numbers: 71.27.+a,71.30.+h,72.10.-d 



I. INTRODUCTION 

The physical nature of scattering processes which con- 
trol transport represents one of the most fundamental 
properties for any material. At the lowest temperatures 
the thermal excitations are few, and elastic impurity scat- 
tering dominates. Raising the temperature introduces 
two basic pathways to modify transport. First, clastic 
scattering can acquire a temperature dependence cither 
through the modified screening of the impurity potential, 
or through dephasing processes**^ This general mech- 
anism encapsulates the physical content of all "quan- 
tum corrections" - both in the diffusive and the ballistic 
regime - predicted within the Fermi liquid framework. 
Indeed, careful and precise experiments have confirmed 
the validity of this physical picture for many good metals 
with weak disorder i Physically, it relies on the existence 
of long-lived quasiparticles within a degenerate electron 
gas. 

The second route comes into play in instances where 
correlation effects due to electron-electron interactions 
are significant. Here, the Fermi liquid regime featur- 
ing degenerate quasiparticles is often restriced to a very 
limited temperature range T <C T* -C Tp, well below 
the "coherence temperature" T*, which itself is much 
smaller than the Fermi temperature Tp. In such materi- 
als, which include rare-earth intermctallicsj 3 -^ many tran- 
sition metal oxides^ and several classes of organic Mott 
systems^ - — a broad intermediate temperature regime 
emerges T ~ T* <g; Tp where inelastic electron-electron 
scattering dominates all transport properties. Such scat- 
tering directly reflects the thermal destruction of Lan- 
dau quasiparticles - a situation describing the demise of 
a coherent Fermi liquid. In these materials, in the rele- 



vant temperature range, the electron-phonon scattering 
is much weaker than the electron-electron one. 

When a material is tuned to the vicinity of any metal- 
insulator transition,— both disorder and electron-electron 
interactions are of a priori importance. But which 
of these two scattering mechanisms - elastic or inelas- 
tic - dominates the experimentally relevant tempera- 
ture range? Answering this question should provide 
important clues as to which of the localization mecha- 
nisms dominate in any given material. Unfortunately, 
experimental systems permitting sufficiently precise tun- 
ing of control parameters are generally rather few. An 
attractive class of systems where a dramatic metal to 
insulator crossover is observed in a narrow parameter 
range is provided by two dimensional electron gases 
(2DEG), such as silicon MOSFETs or GaAs/AlGaAs 
heterostructuresJ^— One of the most striking features 
observed in these systems is the pronounced resistivity 
drop on the metallic side of the transition. While conven- 
tional, relatively weak temperature dependence is found 
at high densities (n ^> n c ), very strong temperature de- 
pendence is found near the critical density n c , roughly 
in the same density range n c < n < 2n c where other 
strong correlation phenomena were observed, e.g. large 
to* enhancement Here, pronounced resistivity max- 
ima are observed at T ~ T max (n), followed by a dra- 
matic resistivity drop at lower temperatures, whose phys- 
ical origin remains a subject of much controversy and 
debate^i 2 . 

In this paper we argue that the electron-electron scat- 
tering dominates the transport in a broad concentra- 
tion and temperature range on the metallic side of 
the metal-insulator transition (MIT) in Si MOSFETS 
and GaAs/AlGaAs heterostructures. This conclusion is 
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reached by: (i) A detailed scaling analysis of the metal- 
lic resistivity curves; (ii) Establishing a similarity in 
the transport properties of the 2DEG and well-studied 
strongly correlated materials near the interaction-driven 
MIT; (iii) Making a comparison of the resistivity curves 
in 2DEG with those in a simple model of the Mott MIT. 
Our conclusions favor the interaction-driven (Wigncr- 
Mott) scenario±2ri£ of the MIT in 2DEG and provide 
a guidance for the development of a microscopical theory 
of incoherent transport in diluted 2DEG. 

The remaining part of the paper is organized as follows. 
Our phenomenological scaling of the experimental data 
is shown in Sec. II, and the analogy with strongly corre- 
lated three dimensional (3D) materials is highlighted in 
Sec. III. The scaling analysis in a simple model of the 
Mott MIT is presented in Sec. IV, and Sec. V contains 
the conclusion and discussion. 



II. SCALING ANALYSIS OF THE RESISTIVITY 
MAXIMA 

The experimental data reveal well defined trends in the 
density dependence of the resistivity maxima, suggesting 
a scaling analysis. While many different scenarios for the 
metal-insulator transition predict some form of scaling, 
its precise features may provide clues to what mechanism 
dominates the transport. 

All of the curves displaying a resistivity maximum 
have an almost identical shape [Fig. Q], strongly sug- 
gesting that unique physical processes are responsible 
for a strong temperature dependence of the resistivity 
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FIG. 2: (Color online) Scaled resistivity as a function of scaled 
temperature for different electron (hole) concentrations, for 
Si MOSFET (a) and GaAs heterostructures (b). The exper- 
imental data are taken from Ref. [U (MOSFETs), Ref. [H 
(p-GaAs/AlGaAs, blue s ym bols), Ref. |24| (n-GaAs/AlGaAs, 
green symbols), and Ref. 12a (p-GaAs, orange symbols). The 
solid line is the scaling function obtained for a simple model 
of the MIT (see Sec. IV). 
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FIG. 1: Resistivity as a function of temperature from the 
experiments on Si MOSFET by Pudalov et al. (Ref. QJ). 



in a large range of concentrations. The resistivity max- 
ima are typically observed at temperatures compara- 
ble to the Fermi temperature, where a physical picture 
of long-lived quasiparticles is no more valid. Comple- 
mentary experimcnt o 11 ' 13 on the same material have re- 
vealed that large effective mass m* enhancements are ob- 
served in the same density range. This behavior is a 
clear signature of strong correlation effects which, in all 
known examples, produce very strong inelastic electron- 
electron scattering in the appropriate temperature range. 
The electron-phonon scattering is negligibly small for 
T < Tp < 10 K3i Since a strongly correlated system 
is typically characterized by a single characteristic en- 
ergy scale T* ~ (m/m*)Tp, we expect the scaling func- 
tion f(x) to assume a universal form, while the scaling 
parameters T max = T* and p max to assume a simple, 
power-law dependence on the effective mass m*. Guided 
by these observations, in this Section we introduce a scal- 
ing ansatz and perform a scaling analysis of the resistivity 
curves in Si MOSFETs and GaAs heterostructures. 
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A. Phenomenological scaling hypothesis 

In accordance to what is typically found in other ex- 
amples of strongly correlated metals with weak to mod- 
erate disorder^ we expect the resistivity to assume an 
additive form, p(T) = p a + Sp(T). Here, p Q is the 
residual resistivity due to impurity scattering, and the 
temperature-dependent contribution Sp(T) is expected 
to be dominated by inelastic electron-electron scattering. 
Based on these general considerations, we propose that 
the temperature-dependent term assumes a scaling form 

Sp(T) = 5p max f(T/T max ), (2.1) 

where Sp max = p max - p a . 



a single dominant transport mechanism. This conclusion 
is even more convincing if we apply the same analysis to 
several different materials, including an ultra high mobil- 
ity GaAs sample, Fig. [UJb). While the diffusive physics 
cannot possible apply in such a broad parameter range, 
we see that the scaling form we propose proves to be 
an extremely robust feature of all available 2D-MIT sys- 
tems. This result is very significant, because disorder 
effects must be significantly weaker in these ultra-clean 
materials, while the interaction effects arc expected to be 
even stronger. 



B. Critical behavior of the Wigner-Mott scaling 
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FIG. 3: (Color online) Tmax normalized to Fermi temperature 

(a) , and maximal resistivity 8 pmax = pmax — Po in units nh/e 2 

(b) , as a function of a reduced density. The data are taken 
from Refs. HHHM 

To test this phenomenological scaling hypothesis, we 
perform a corresponding analysis of experimental data 
in several systems displaying 2D-MIT. We start with the 
Si MOSFET dataiS analyzed in Ref. [H We concen- 
trate on metallic curves below the separatrix C. In the 
range of concentrations, 0.83 < n < 1.10, the resistiv- 
ity curves have a clear maximum, and nicely collapse 
with the proposed scaling ansatz, Fig. [2ja). In fact, we 
can use the scaling ansatz to collapse also the data for 
1.21 < n < 1.75, where T max and p max are determined 
from the least square fit to the scaling curve. Clearly all 
eight resistivity curves belong to the same family (have 
the same functional form) , and thus must be explained by 
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FIG. 4: (Color online) (a) Maximum resistivity Sp, n ax = 
pmax - po as a function of T ma x- (b) T,„ax as a function 
of inverse effective mass m* . m is the band mass in Si MOS- 
FETs. The data are taken from Refs. [2l| and [H 

Having demonstrated data collapse, we are now in a 
position to examine the critical behavior of the rele- 
vant crossover scale. We thus examine the behavior of 
T max and p max as a function of reduced concentration 
(n — n c )/n c and effective mass to* (as determined by 
complementary experiments). 

For different realizations of 2DEG, T max shows approx- 
imately power law dependence on the reduced concentra- 
tion [Fig.[3ja)] and even the exponents are similar. T max 
in our physical picture has a clear physical interpretation 
as a coherence temperature - the temperature when the 
inelastic electron-electron scattering time becomes com- 
parable to h/Ep, leading to incoherent transport. The 
resistivity maximum, however, shows less universal form. 
It varies a lot in different physical systems. This does not 
come as a surprise since the resistivity shows nonuniver- 



4 



sal features also in three dimensional strongly correlated 
materials near the Mott transition. We discuss in detail 
the analogy with the Mott systems in Sees. Ill and IV. 

In a Si MOSFET the resistivity maximum 8p max = 
Pmax—Po shows power law dependence on T max in a fairly 
broad concentration range [Fig. EJa)]. We further ana- 
lyze the critical behavior for Si MOSFET using the data 
for the effective mass as determined by Shashkin et alJ£ 
from magnetorcsistance measurements in a parallel mag- 
netic field. We find that T max is inversely proportional 
to the effective mass m*. This behavior is typical to all 
systems near the Mott MIT, where the coherence tem- 
perature is inversely proportional to the effective mass, 
as a landmark of strong correlations. 

C. Breakdown of the diffusion mode scaling 



clastic scattering rate H/t. According to this picture, a 
different (ballistic) mechanism for transport is expected 
outside the diffusive regime, presumably leading to a dif- 
ferent temperature dependence, so the proposed scaling 
no longer holds. This analysis was applied to the experi- 
mental data of Refill, but was accordingly restricted to 
only three densities closest to the transition. Indeed, if 
the scaling formula is applied in a broader range of con- 
centrations, the resistivity curves clearly do not collapse 
[Fig. [5]. While the Fermi liquid renormalization group 
calculations are very important in order to answer a fun- 
damental question of necessary conditions for a true MIT 
at zero temperature, our analysis emphasizes that the un- 
derstanding of various diluted 2DEG in a broad range of 
parameters requires the physics beyond the conventional 
Fermi liquid framework. 



We have successfully collapsed resistivity curves in a 
broad temperature and concentration range and for sev- 
eral physical realizations of 2DEG. The physical pic- 
ture behind the proposed scaling is that the 2D MIT 
is an interaction-driven (Wigner-Mott) MIT^~— and 
that the dominant temperature dependence in the re- 
sistivity originates from strong electron-electron scatter- 
ing. Another proposed scenario envisions disorder as the 
principal driving force for localization t 21 i 22 while the in- 
teractions are most important above the critical density 
and at low temperatures, where they suppress the ten- 
dency to localization. An appropriate theory, based on a 
Fermi liquid framework^ has predicted that a resistivity 
maximum should be observed on the metallic side, with 
the resistivity assuming the scaling form 



P(T)/Pr> 



f[Pmax ^(T/^max)] 



(2.2) 



Here f(x) is a universal scaling function predicted by 
theory. The authors point out, though, that this pre- 
diction is expected to be valid only within the diffusive 
regime, where the thermal energy ks T is smaller than the 
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III. SCALING IN 3D MATERIALS 

The strong temperature dependence of resistivity is a 
well known feature of many strongly correlated materials. 
A pronounced resistivity maximum is observed in heavy 
fcrmionsj^ and charge-transfer organic salts where 
the correlation strength is tuned by applying an exter- 
nal pressure. The essential mechanism of transport in 
these materials relies on strong inelastic electron-electron 
scattering, and the Fermi liquid behavior is restricted to 
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FIG. 5: (Color online) Resistivity as a function of temperature 
scaled as in Ref.— . Red solid line is the calculated scaling 



curve. 



FIG. 6: (Color online) Scaled resistivity curves for UBei3 
(a) and k-(ET)2Cu2(CN)3 (b), for different external pressure. 
The data are taken from Refs. and 0- 
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the lowest temperatures. As the temperature increases, 
the electron mean free path becomes comparable to, or 
smaller than the lattice spacing, and the transport be- 
comes incoherent. The electron-phonon scattering is here 
much weaker than the electron-electron one. The tem- 
perature of the resistivity maximum can be taken as a 
definition of the coherence temperature T* . It is inversely 
proportional to the effective mass, and much smaller than 
the bare Fermi temperature, T* ~ (mb/m*)Tp. The 
same scaling ansatz as given by Eq. (|2.1[) was used to 
collapse the resistivity curves for CeCug already in an 
early paper by Thompson and Fiskj^ 

Here we illustrate the similarity in transport proper- 
ties of these systems and 2DEG by scaling the resistivity 
data for heavy fermion UBei3 from Ref. [Fig. [6(a)], 
and for a charge-transfer conductor k-(ET)2Cu2(CN)3 
[Fig. E(b)]. The collapse of the resistivity curves is ex- 
cellent for UBei3, and well-defined trends are seen in 
K — (ET) 2 Cu 2 (CN)3. Remarkable similarity in resis- 
tivity curves in such diverse physical systems like Si 
MOSFETs, GaAS heterostructures, heavy fermions and 
charge-transfer organic conductors is, in our view, a man- 
ifestation of the same physical processes in the vicinity 
of the interaction-driven MIT. 



IV. SCALING IN THE MICROSCOPIC MODEL 
OF THE INTERACTION-DRIVEN MIT 

Having phenomenologically established precise and 
well defined scaling behavior of the experimental curves 
on the metallic side of the 2D MIT for temperatures near 
T* , we now address its microscopic origin. More pre- 
cisely, we would like to understand just how robust this 
result is. Does it depend on subtle details describing the 
interplay of disorder and interactions of 2DEG materials, 
as suggested in Ref. [27], or is it a generic feature of strong 
correlation near interaction-driven MIT. To answer this 
important question we deliberately focus on the sim- 
plest microscopic model for interaction-driven MIT: The 
clean single-band Hubbard model at half-filling. Accu- 
rate and quantitatively precise results can be obtained for 
temperature-dependent transport for this model within 
the DMFT approximation^ While the DMFT repro- 
duces Fermi liquid behavior at the lowest temperatures, 
it is particularly useful in the studies of " high tempera- 
ture" incoherent transport. Results of such calculation, 
obtained by the Continuous Time Quantum Monte Carlo 
(CTQMC) impurity solve r 29 i 30 followed by the analytical 
continuation by the Maximum Entropy Method 3 -^, can 
be analyzed using precisely the same scaling procedure 
we proposed for experimental data. We concentrate on 
the metallic phase of the Hubbard model with the inter- 
action parameter U smaller than the value at the criti- 
cal end-point U c . The resistivity curves [Fig. [7(a)] have 
qualitatively the same form as in 2DEG. The resistivity 
sharply increases with temperature, reaches a maximum 
and than decreases. The temperature of resistivity max- 
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FIG. 7: (Color online) (a) Resistivity as a function of temper- 
ature for several interaction strengths in the half-filled Hub- 
bard model solved within the DMFT. The resistivity is nor- 
malized to the Mott limit value, which corresponds to the 
scattering length of one lattice spacing, (b) Scaled resistivity 
curves. 



imum decreases as the system approaches the MIT. 

Most remarkably, precisely the same scaling form as in 
2DEG is found to describe all resistivity curves close to 
the Mott transition [Fig. [7f b)] . In addition, we find that 
the scaling parameters T max and p ma x again display a 
power law dependence on the effective mass [Fig. [5] , and 
even the exponents are similar. Finally, we contrast the 
DMFT scaling function with that obtained from 2DEG 
experiments. We find surprisingly accurate agreement 
between the DMFT prediction for the scaling function 
f(x) and experimental data on all available materials 
[Fig. [2]. We emphasize, however, that our scaling hy- 
pothesis is valid only in the metallic phase for U < U c 
and for temperatures comparable to T* ~ 1/m*. It 
should be contrasted with the scaling near the critical 
end-point (U c , T e ) ) 32 i 33 or the proposed quantum critical 
scaling in the high-temperature regime above the critical 
end-point i2i 

We should point out that for this model, the proposed 
resistivity scaling is not valid at the lowest temperatures 
T <C T max , deep within the Fermi liquid region: Ac- 
cording to the Kadowaki- Woods relation, here p « AT 2 
where A ~ 1/m* 2 ~ 1/T^ ax , and the scaling is violated 
if the resistivity is scaled by p ma x- For T > 0.3T max the 
collapse of the resistivity curves is excellent, [Fig. [7(b)], 
and we define the DMFT scaling curve for this temper- 
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FIG. 8: (Color online) (a) Maximum resistivity as a function 
of the corresponding temperature from the DMFT solution 
of the Hubbard model, (b) T max as a function of the inverse 
effective mass. 



ature range. This is also the reason for the deviations 
in the scaling in Fig. [Bfb) for K-organics, the materials 
whose properties are described remarkably well within 
the Hubbard modeli^ In the Anderson lattice model, on 
the other hand, the resistivity maximum does not change 
much near the MIT and it saturates approximately to 
the value which corresponds to the scattering length of 
one lattice spacing (Mott limit). In this case our scal- 
ing ansatz is valid in the whole temperature range up to 
T = 0, 36 and the collapse of the resistivity curves seen in 
the experiments is excellent [Fig. [5[a)]. 

Microscopic theory of the 2DEG should also in- 
clude nonlocal correlations which are neglected in 
a simple DMFT approach. A more realistic ex- 
tended Hubbard model displays a two-stage Wigner-Mott 
localization! 17 ' 18 The metal-insulator transition in this 
model is found in the region with already developed 
nonlocal charge correlations. In the immediate critical 
regime, the critical behavior can be represented by an 
effective Hubbard model, partially justifying the success 
of the present modeling. The existence of a coherence 
scale T* which vanishes at the onset of charge order is 
also found in the 2D extended Hubbard model solved by 
finite- T Lanczos diagonalization.— This result is relevant 
for quarter-filled layered organic materials, which further 
supports the importance and generality of the ideas pre- 
sented here. 



V. CONCLUSION AND DISCUSSION 

In this paper we argued that the emergence of resis- 
tivity maxima upon thermal destruction of heavy Fermi 
liquids should be regarded as a generic phenomenon in 
strongly correlated systems. We demonstrated that the 
resulting family of resistivity curves typically obeys a 
simple phenomenology displaying scaling behavior. Our 
detailed model calculations show that all of the qualita- 
tive and even quantitative features of this scaling phe- 
nomenology are obtained from a microscopic model of 
heavy electrons close to the Mott metal-insulator transi- 
tion. We should stress, however, that the proposed scal- 
ing behavior obtains - both in our theory and in exper- 
iments - only within the metallic regime not too close 
to the transition and the temperature regime around the 
resistivity maxima. In contrast, earlier experiments fo- 
cused on the immediate vicinity of the metal-insulator 
transition, where different " quantum critical" scaling was 
foundi 10 ' 37 ~ — Remarkably, precisely such behavior was 
also found in very recent studies of quantum critical 
transport near interaction-driven transitions^ but this 
was identified in a different parameter regime than the 
one studied in the present paper. 

Our results provide compelling evidence that sev- 
eral puzzling aspects of transport in low density two- 
dimensional electron gases in zero magnetic fields can 
be understood and explained within the Wigner-Mott 
scenario of strong correlation*^— This physical picture 
views the strong correlation effects in the low density 
2DEG as the primary driving force behind the transi- 
tion, and additional disorder effects as less significant, 
secondary processes. In the Wigner-Mott picture the in- 
sulator essentially consists of interaction-localized mag- 
netic moments. Remarkably, the magneto-capacitance 
measurements of Prus et al^ show that the behavior 
characteristic of localized magnetic moments, x(T)/n 
g/jig/T, is seen near the critical density, while only 
weak Pauli-likc temperature dependence was observed at 
higher density. Very recent experiments on Si MOSFETs 
find that the thermopower diverges near the MIT4£ The 
authors argue that divergence of the thermopower is not 
related to the degree of disorder and reflects the diver- 
gence of the effective mass at a disorder-independent 
density, behavior that is typical in the vicinity of an 
interaction-induced phase transition. Additional hints 
supporting this physical picture of 2D MIT are provided 
by existing first principle Quantum (diffusion) Monte 
Carlo results for the low density 2DEG of Ceperley^ 1 - and 
othersi 27 ' 42 ' 43 These calculations find that the correlated 
metallic state has an "almost crystalline" structure, thus 
having a very strong short range charge-order (as seen, 
for example, in the density correlation function). 

Within the physical picture that we propose, 
the inelastic electron-electron scattering takes central 
stage ; 44 ' 45 in contrast to disorder-dominated scenar- 
ios, where the interaction effects mainly introduce 
the temperature dependence of elastic electron-impurity 
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scattering^ The two physical pictures describe two com- 
pletely different scattering processes, which are expected 
to be of relevance in complementary but in essentially 
non-overlapping parameter regimes. Indeed, inelastic 
scattering dominates only outside the coherent Fermi- 
liquid regime, which in good metals happens only at fairly 
high temperatures. In strongly correlated regimes that 
we consider, the situation is different. Here the Fermi 
liquid coherence is found only at very low temperatures 
T < T* <C Tp, behavior which is generally observed in 
all system with appreciable effective mass enhancement. 
The results presented in this paper provide precise and 
detailed characterization of this incoherent regime, re- 
vealing remarkable coincidence of trends observed in the 
experiment to those found from the Wigner-Mott picture 
of the interaction-driven metal-insulator transition. Our 
scaling ansatz is proposed based on the physical argu- 
ments and the experimental data. While consistent with 
simple model calculations for strongly correlated elec- 
tronic systems, our work does not directly address spe- 
cific microscopic mechanism responsible for current dissi- 



pation, a process that in 2DEG systems should be facili- 
tated by impurities and imperfections 4£ Still, it provides 
very strong motivation to develop a more realistic mi- 
croscopic theory of incoherent transport in the strongly 
correlated regime of diluted 2DEG. This important task 
remains a challenge for future work. 
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